Elastocaloric effect associated with the martensitic transition in shape-memory alloys 
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The elastocaloric effect in the vicinity of the martensitic transition of a Cu-Zn-Al has been studied by induc- 
ing the transition by strain or stress measurements. While transition trajectories show significant differences, 
the entropy change associated with the whole transformation (AS t ) is coincident in both kind of experiments 
since entropy production is small compared to A St. The values agree with estimations based on the Clausius- 
Clapeyron equation. The possibility of using these materials for mechanical refrigeration is also discussed. 

PACS numbers: 65.40.Gr, 75.30.Sg, 81.30.Kf 



Caloric effects are expected to occur under the application 
of an external field to a given material. The elastocaloric ef- 
fect HI] is the mechanical analogue of the magnetocaloric ef- 
fect that has received considerable attention in the recent years 
owing to its potential use for environmentally friendly refrig- 
eration 10] ■ The magnetocaloric effect is related to the isother- 
mal change of entropy or the adiabatic change of tempera- 
ture that takes place within a material when a magnetic field 
is applied or removed. This effect originates from the cou- 
pling between the magnetic sublattice and an externally ap- 
plied magnetic field and thus occurs in any magnetic material. 
A large effect is expected in the vicinity of field-induced, first- 
order phase transitions where large entropy changes should 
occur 01, By analogy, the elastocaloric effect is defined as 
the isothermal change of entropy or the adiabatic change of 
temperature that takes place when a mechanical field (stress) 
is applied or released in a given material. Indeed, this effect 
is expected to be a consequence of the coupling between an 
external applied stress and the lattice. Continuing with the 
analogy, a large elastocaloric effect is also foreseen in sys- 
tems undergoing stress-induced, first-order phase transitions. 
Good candidates to show this effect are shape-memory alloys. 
These materials undergo a diffusionless purely structural tran- 
sition from a cubic to a lower symmetry phase that can be 
stress induced 0]. Actually, shape-memory properties are re- 
lated to this transition and refer to the ability of these systems 
to remember their original shape after severe deformation ^ . 

In contrast to magnetism, instead of controlling the ap- 
plied stress (or force) which is the variable thermodynami- 
cally equivalent to the magnetic field, in mechanical experi- 
ments, the system is usually driven by controlling the strain 
(generalized displacement) which is the conjugated variable 
to the stress in the way that magnetization is the conjugated 
variable of the magnetic field. In magnetic systems, due to the 
difficulty in controlling magnetization, magnetization-driven 
experiments aimed at studying the magnetocaloric effect have 
not, to our knowledge, been reported. Thus, comparing results 
from both field/stress-driven and magnetization/strain-driven 
experiments is of general interest since constraining the (gen- 
eralized) displacement prevents free motion of the interfaces 
and therefore field/stress fluctuation can occur. Indeed, this 
is especially relevant in systems undergoing a phase transi- 



tion leading to macroscopic instability. The study of mechan- 
ical systems naturally opens up the possibility of performing 
such kind of experiments. In a recent work ]6[] we have shown 
that metastable trajectories exhibit a strong dependence on the 
driving mechanism. In particular, strain-driven trajectories are 
characterized by the occurrence of re-entrant behaviour and 
lower dissipation than stress-driven trajectories. 

The present letter is aimed at studying the elastocaloric 
effect in the vicinity of the martensitic transition in a 
Cu68.13Zn15.74Ali6.13 single crystal (molar volume = 7.52 
cm 3 mol _1 ) in both the stress-driven and the strain-driven 
modes. The sample was mechanically machined with cylin- 
drical heads and the body has flat faces 35 mm long, 4 mm 
wide and 1.4 mm thick. Its axis is close to the [100] crystal- 
lographic direction of the cubic phase. This orientation pro- 
vides high transformation strains for relatively low uniaxial 
applied tensile loads along this direction thus preventing the 
occurrence of irreversible plastic effects. The sample was con- 
veniently heat treated to ensure that it was free from internal 
stresses and that the order state and vacancy concentration are 
close to their ground state values JjJ. In the absence of any 
applied stress, the crystal undergoes a martensitic transition 
from a cubic L2i phase to a multivariant 18R martensite at 
T M = 234 K. 

For strain-driven experiments we used an Instron 4302 
screw driven tensile machine in which the elongation is the 
control parameter. For stress-driven experiments we used a 
machine which enables control of the force applied to the sam- 
ple while elongation was continuously monitored [6]. The 
machine applies a dead load to the sample which can be in- 
creased or decreased at a well controlled rate. A cryofurnace 
(with temperature stability ±0.1 K) can be adapted to both 
devices. All experiments have been performed at low rates [~ 
0.3 mm min" 1 (strain-driven) and ~ 5 MPa min -1 (stress- 
driven)]. 

Let us now consider a thermodynamic system described by 
variables {X, Y, T}, where X is a generalized displacement 
and Y is the corresponding conjugated field (X and Y have 
the same tensorial order), and T is the temperature. A change 
in the generalized displacement gives rise to a caloric effect. 
If this change is induced by an isothermal change AY of its 
conjugated field, the caloric effect must be quantified by the 
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corresponding induced entropy change which is given by 

where the generalized Maxwell relation (9S/<9Y)t = 
(<9X/<9T)y has been taken into account. When the field is 
an uniaxial tensile stress a for which the corresponding strain 
(or relative elongation along the direction of the applied force) 
is e, the induced isothermal entropy change defining the elas- 
tocaloric effect, is given by 



AS(0 



(2) 



This expression is formally analogous to the expression giv- 
ing the field-induced entropy change which defines the mag- 
netocaloric effect in a magnetic system J2]. If instead of a the 
controlled variable is e, the entropy change corresponding to 
an isothermal variation of the strain from to e is given by 
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where the Maxwell relation (dS/de)r = —(da/dT) e has 
been used in this case. Of course, if e is the strain corre- 
sponding to the stress a, in strict equilibrium AS(0 — > a) — 
AS(0 -> e). 

Now assume a system subjected to an applied uniaxial 
stress a that undergoes a first-order structural (martensitic) 
phase transition in equilibrium at a temperature T t . The transi- 
tion is in this case characterized by discontinuities in variables 
such as strain and entropy that are thermodynamically conju- 
gated to the intensive variables stress and temperature. In the 
vicinity of the transition the following behaviour of the strain 
is thus expected, e(T,cr) = e + Ae T[(T t {cr) - T)/AT], 
where T is a shape-function and AT is a measure of the tem- 
perature range over which the transition spreads. In strict 
equilibrium, AT — > so that T approaches the Heaviside 
step function. Using expression (f2]i and assuming that e and 
Ae are constant, in this equilibrium case the elastocaloric ef- 
fect in the vicinity of the transition is given by 



AS(0 -> a) 







forTG [T t (0),T t (cr)] 
forT£ [T t {0),T t (a)} 



(4) 



where a = dTt/da is assumed to be constant. Taking into 
account the Clausius-Clapeyron equation, a = —Ae/ASt, 
where A St is the transition entropy change. Therefore, as 
expected, AS(0 -> a) = AS t , and AT = T a - T t (0) = cut. 
Indeed, the same result is obtained in this case from Eq. 0. 

Actually, these transitions are spread over a small range 
of a, and, more importantly, they display hysteresis, which 
reflects the existence of non-equilibrium dissipative effects. 
In this case, taking into account the Clausius inequality, 
§ Sq/T < 0, for an isothermal process, the entropy change 
must satisfy, AS(0 — > a or e) = ^ + Si, where Si > is the 



entropy production and is expected to depend on the actual tra- 
jectory followed by the system (strain- or stress-driven in our 
case). The preceding equation indicates that an estimation of 
the entropy change based on expression (0 should differ from 
an estimation based on direct calorimetric measurements. For 
the system of interest here, as hysteresis is quite small, estima- 
tions of the entropy change based on stress-driven and strain- 
driven curves should provide reasonably good estimations of 
the elastocaloric effect. 

In Fig. [T]we show stress-strain curves obtained in the stress- 
driven case at selected temperatures (well above the transition 
temperature at zero-stress) across the martensitic transition in 
the studied Cu-Zn-Al single crystal. The shift of the transition 
to higher stresses with increasing temperature is clearly seen. 
Comparison of curves corresponding to loading and unloading 
shows that the transition occurs with weak hysteresis of about 
10 MPa. From these curves the elastocaloric effect (stress- 
induced entropy change) has been obtained by numerically 
computing the integral in Eq. ©. The stress-induced entropy 
change (elastocaloric effect) is shown in Fig. [2] 

It is interesting to point out that the maximum stress- 
induced entropy change (which corresponds to the whole en- 
tropy change of the martensitic transformation) remains al- 
most constant (thus independent of temperature and applied 
stress) over a very broad temperature range [7]. This is usu- 
ally not the case for the magnetocaloric effect in field-induced, 
first-order phase transitions for which large AS(T) is only 
obtained in a relatively narrow temperature interval (which 
depends on the applied field). Such a difference is a conse- 
quence of the fact that here tensile experiments are performed 
at temperatures well above (> 60 K) the transition temper- 
ature at zero stress. This is possible due to the strong de- 
pendence of the transition stress with temperature. The up- 
per bound is imposed by the elastic limit of the cubic phase. 
By contrast, in magnetic experiments the field is always ap- 
plied close to the transition temperature at zero field. Indeed, 
at higher temperatures metamagnetic transitions are difficult 
to be induced either because intense magnetic fields that are 
too large are required or simply because this leads the system 
above the critical point, where no transition occurs. The con- 
tinuous lines in Fig. |2]are fits which assume a shape-function 
T{x) = tanh - (a;). In the strain-driven case we can proceed 
similarly starting with strain-stress curves recorded at selected 
temperatures. Fig. [3]displays the strain-stress curves obtained 
at selected temperatures. The strain-induced entropy change 
(elastocaloric effect) as a function of T has then been obtained 
using Eq. ([3). Results are shown in Fig. [4] The variation 
of AS with temperature reflects a small variation of Ae with 
temperature. 

From the previous results the estimated entropy change 
corresponding to the whole transition is ASt = —1-20 ± 0.15 
J/mol K consistent in both strain- and stress-driven cases. It 
is interesting to compare this result with estimations based on 
the use of the Clausius-Clapeyron equation. From the stress- 
driven curves in Fig. Q~|we can define a transition stress, a t , 
at each temperature, as the stress at the inflection point of 
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FIG. 1: (Color online) Stress-strain curves at selected temperatures. 
Continuous curves correspond to loading runs (from right to left) at 
T = 310.5, 309.4, 307.5, 306.0, 305.4, 303.1, 302.0, 299.9, 298.8, 
297.9, 296.1, 295.3 and 294.6 K. The left discontinuous curve cor- 
responds to the unloading branch at T — 294.6 K. It illustrates an 
example of stress-induced hysteresis loops. The inset shows the tran- 
sition stress as a function of temperature. The line is a linear fit 
a = 2.01(T - 242). 
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FIG. 3: (Color online) Strain-stress curves obtained at selected tem- 
peratures. Continuous curves correspond to loading runs at (from top 
to bottom): T = 307.8, 303. 1, 297.8, and 295.0 K The lower discon- 
tinuous curve corresponds to the unloading branch at T = 295.0 K. 
It illustrates an example of a strain-induced hysteresis loop. 
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FIG. 4: (Color online) Strain-induced entropy change (elastocaloric 
effect) at selected values of e from 0.01 to 0.08. 



FIG. 2: (Color online) Stress-induced entropy change (elastocaloric 
effect) at selected values of a ranging from 105 to 143 MPa. The 
continuous lines are fits based on a the model T[x) = tanh - (x). 



the isotherms. The strain change, Ae, at the transition can 
also be estimated from these curves. The transition entropy 
change is then obtained as ASt — (dat/dT)As. In the in- 
set of Fig. Q]we show a t vs. T. We see that a t rises lin- 
early with increasing T, with a slope dat/dT = 2.01 MPa/K. 
Taking an average value Ae = 0.080 ± 0.005 we obtain 
AS t = -1.21 ± 0.05 J/mol K. This value is compatible with 
the estimations based on elastocaloric effect. It is interest- 
ing to compare this value with a calorimetric measurement. 
To this end, we have carried out calorimetric measurements 
of the transition entropy change at zero stress, using a small 
specimen cut from the same original ingot. These experiments 
give a value ASt = —1.37 ± 0.10 J/mol K, slightly higher 
(as an absolute value) than that derived from mechanical ex- 
periments. Such a difference should not, in principle, be at- 
tributed to non-equilibrium effects which have been estimated 
from the area of the hysteresis loops to be of the order of 0.01 
J/mol K in the stress-induced case and less than 0.001 J/mol 
K in strain-induced experiments. However, it must be pointed 
out that calorimetric measurements are performed in the ab- 



sence of an external field and thus a multivariant martensite is 
reached. Kinematic constraints which occur during the trans- 
formation and yield extra dissipative effects could explain the 
difference. This is consistent with the fact that the transition 
temperature in the multivariant case (234 K) is significantly 
lower than the extrapolation to zero stress of the a vs. T curve 
(242 K) which corresponds to the transition to a single variant 
martensite (see inset of Fig. [T). 

The elastocaloric effect has been reported in the Fe-Rh al- 
loy JmI ■ In this case the entropy change has a magnetic origin 
associated with a re-orientation of spins taking place at the 
first-order antiferromagnetic to ferromagnetic magnetostruc- 
tural transition which is induced by the application of a tensile 
stress due to strong magnetoelastic coupling. Similar effects 
(termed the barocaloric effect) have been reported for a num- 
ber of rare-earth compounds subjected to uniaxial pressure 
Jioll . In the present study, however, the elastocaloric effect is 
associated with a purely structural transition. In this case, the 
entropy change is predominantly vibrational and originates 
from the very low energy TA2 transverse ([110] propagation 
and [1 10] polarization) phonon branch of the cubic phase @]. 
It is also worth noticing that in Fe-Rh the entropy change is 
positive when the stress is isothermally applied, and thus the 
sample cools down when the stress is adiabatically applied. 
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This corresponds to an inverse elastocaloric effect which is 
the analog of the inverse magnetocaloric effect reported in 
Heusler martensitic alloys ll ill . Notice that inverse effects are 
only possible when there is strong coupling between magnetic 
and structural degrees of freedom. 

The adiabatic temperature change associated with the elas- 
tocaloric effect can be estimated as AT ~ ~^AS, where C 
is the specific heat which is assumed to be stress independent. 
For Cu-Zn-Al, the value for the specific heat close to room 
temperature is approximately 25 J/K mol in both martensitic 
and cubic phases IU2I1 . For an adiabatic stress drop of about 
10 MPa (for instance from 1 10 MPa to 100 MPa at 300K, see 
Fig. [TJ the maximum expected temperature change is 15K. 
Notice that this value is orders of magnitude larger than the 
typical values in elastic solids far from any phase transition. 
The effect of strain rate on the martensitic transition of several 
Cu-based shape-memory alloys was studied by means of di- 
rect measurements of the temperature changes associated with 
the transition. For Cu-Al-Ni HI] and Cu-Zn-Sn J3 crystals 
(with the same martensitic structure as our Cu-Zn-Al crystal), 
the measured temperature changes at high strain rates (close to 
the adiabatic limit) are 14 K and 12 K, respectively, for a strain 
rate of 25 rnin -1 . These temperature values compare well to 
the value indirectly computed here for Cu-Zn-Al. On the other 
hand, present values for Cu-Zn-Al are comparable to those re- 
ported for other elastocaloric materials undergoing first-order 
magnetostructural phase transitions: maximum changes of 8.7 
K and 14 K are computed for Fe-Rh H and Eu-Ni-Si-Ge ifioll , 
respectively. 

The elastocaloric effect associated with the martensitic 
transition in a Cu-Zn-Al single crystal has been studied. It is 
formally equivalent to the magnetocaloric effect in a magnetic 
system and in our case it should be comparable to the magne- 
tocaloric effect in the vicinity of a first-order metamagnetic 
transition. However, while the magnetocaloric effect is al- 
ways determined by the field inducing the metamagnetic tran- 
sition, for the mechanical case it has been possible to obtain 
the isothermal entropy change by inducing the structural tran- 
sition using both strain and stress. While the transition path is 
essentially different in both cases, the corresponding isother- 



mal entropy changes are the same to within errors. This is due 
to the fact that hysteresis is small independent of the driving 
mechanism. For practical applications of caloric effects, the 
refrigerant capacity is a central parameter to be considered ■ 
It is defined as K = J AT AS(T)dT ~ A SAT = -AeAa, 
where Aa is the change of a necessary to change the transi- 
tion temperature by AT (AS" and Ae are assumed constants). 
The interest in our case is that Act can be chosen in a broad 
range of values which opens up interesting opportunities in 
refrigeration applications based on the elastocaloric effect. 
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